Nanopatterning of Si surfaces by normal incident 
ion erosion: influence of metal incorporation on 
surface morphology evolution 
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Abstract. The morphology evolution of Si (100) surfaces under 1200 eV Ar + 
ions bombardment at normal incidence with and without metal incorporation is 
presented. The formation of nanodot patterns is observed only when the stationary Fe 
concentration in the surface is above 8 x 10 14 cm~ 2 . This is interpreted in terms of an 
additional surface instability due to non-uniform sputtering in connection with metal 
enrichment at the nanodots. At low metal concentration, smoothing dominates and 
pattern formation is thus inhibited. The transition from a fc~ 2 to a fc~ 4 behavior in 
the asymptotic power spectral density function supports the conclusion that ballistic 
smoothing and ion-enhanced viscous flow are the two dominant mechanisms of surface 
relaxation. 
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1. Introduction 

Low energy ion sputtering is a promising candidate for the fabrication of self-ordered 
nanostructures on both crystalline and amorphous surfaces [Tj. It was observed 
experimentally that ion sputtering at normal incidence or at oblique incidence can 
generate correlated dot or ripple patterns on metals [2], semiconductors [31 EJ [5] and 
insulators [6]. The formation mechanism has been described theoretically by the 
competition between morphology induced surface instability and surface relaxation 
processes pj El El [TO] . 

Recently, it was reported by Ozaydin et al [H] that molybdenum incorporation 
triggers the formation of ordered dot arrays on Si surfaces under the bombardment 
with Ar + ions at normal incidence while no correlated structures are generated in the 
absence of the impurities. Furthermore, Hofsass et al [12] found that co-deposition of 
Au, Ag and Pt surfactants generates novel patterns and nanostructures on Si surfaces 
during ion sputtering which are absent without co-deposition. Moreover, Sanchez- 
Garcia et al [13] observed a transition from hole to dot patterns on Si by tuning the 
amount of metal incorporation during ion sputtering at normal incidence. Additionally, 
Shenoy et al [H] proposed that preferential sputtering can generate modulation both in 
surface morphology and composition. These findings add a new aspect to the present 
understanding of the development of surface topography during ion sputtering and 
suggest more detailed studies which address the involved roughening and smoothing 
mechanisms. Especially quantitative data on the influence of the metal incorporation 
on the morphology evolution are still missing. 

In this contribution, a comparative experiment is described, which proves the 
important role of metal incorporation in triggering dot pattern formation on Si(100) 
surfaces and allows us to determine the Fe concentration at which the morphology 
evolution changes from roughening to smoothing.. Furthermore, the morphology 
evolution in the absence of metal incorporation clearly indicates that the relaxation 
process includes both k 2 - and /c 4 -dependent smoothing mechanisms, with k being the 
surface wavenumber. 

2. Experiment 

The experiments were performed in an ultra-high- vacuum chamber with a base pressure 
of 1 x 10~ 8 mbar. For generation of a low-pressure plasma, an inductively coupled RF 
plasma (ICP) generator is attached using a planar coil behind a quartz window of 150 
mm diameter, which is immersed into the chamber [15]. At an Ar pressure of 4xl0 -3 
mbar, the absorbed RF power was 100 W at a frequency of 13.56 MHz. Ions were 
extracted from the plasma by applying a 1200 eV negative DC bias to a flat sample 
holder which is mounted opposite to the quartz window at a distance of 10 cm. The 
electric field geometry and the low pressure guarantee an almost ideal normal incidence 
of the ions [16]. The resulting current density was 200 fiA cm -2 at the sample. During 
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Table 1. Surface contents of Ar, Fe, and Cr as obtained from RBS measurements. 







Contents (10 15 


cm 2 ) 


Sample 


Mount 


Ar Fe 


Cr 


A 


Stainless Steel 


1.63±3% 1.53±3% 


0.40±13% 


B 


Si coverage 


1.64±3% 0.10±50% 


0.01±100% 



operation, the sample temperature remained below 100 °C. It is well established that 
the ion induced pattern formation on Si does not depend significantly on temperature 
between -20 °C and -150 °C [13]. 

Commercially available epi-polished Si(100) wafers (n-type, ~10 fi-cm) were used as 
samples. The metal incorporation was switched on and off by using two different masks 
that also served to fix the samples on the sample holder. Metal atoms sputtered from a 
stainless steel mask may be re-deposited on the sample surface due to partial ionization 
in the plasma. Alternatively, the steel mask was covered by an evaporated layer of Si 
of ~1 fim thickness, so that the metal incorporation is largely suppressed. For post- 
sputtering experiments, the Si(100) surfaces were pre-patterned by ion bombardment 
with the same parameters as stated above, removed from the chamber and re-mounted 
with the Si covered mask for subsequent irradiation at different fluences. 

The surface topography was analyzed ex situ by atomic force microscopy (AFM) in 
tapping mode. The surface metal content was analyzed by Rutherford backscattering 
spectrometry (RBS) using a beam of 1.7 MeV He + ions at normal incidence. Ions 
scattered at an angle of 170° were analyzed using a Si surface barrier detector. 

3. Results and discussion 

3.1. Influence of metal incorporation on dot pattern formation 

Different samples, which will be denoted as A and B in the following, were fixed with a 
bare stainless steel mask and the Si covered mask, respectively, as described above. They 
were irradiated during 15 minutes each, corresponding to an ion fluence of 1.1 x 10 18 
cm . 

Table CD shows the contents of three representative elements on both samples after 
irradiation. The Ar content is not influenced by the amount of the metal impurities 
while the amount of Fe and Cr is reduced on sample B by more than one order of 
magnitude. The relative amounts of Fe and Cr correspond to the atomic fraction of Cr 
in SS 304 (~18 %). Figure Q] shows the power spectral density (PSD) functions obtained 
from the AFM images (inset) of the two samples. An ordered nanodot array is formed 
only in the condition of high metal incorporation, with a mean roughness of 0.85 nm. 
The presence of the peak in the PSD function reveals a characteristic periodicity of 50 
nm. In contrast, no correlated structures are observed on sample B at a roughness of 
0.23 nm. 
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Figure 1. AFM images (inset) and the corresponding PSD functions from Si(100) 
samples after sputtering in the presence and absence of metal incorporation (A and B, 
respectively). Two lines with slopes -2 and -4 are plotted for comparison (see text). 



Basically, the surface morphology evolution results from an interplay between 
surface instabilities and relaxation processes. According to Bradley and Harper (BH) 
[7], the morphology evolution can be described by a partial differential equation of the 
surface height /i(x, t) where x denotes the two-dimensional lateral surface coordinate 
and t the time. At normal incidence, the laterally isotropic BH equation reads, 

^<M = _ Wo + v ^ h _ dv 4 ^. (1) 
at 

Vq is the erosion rate of the unperturbed planar surface, v and D are parameters 
which depend on the experimental conditions. The second term on the right hand side 
represents the curvature dependent roughening and smoothing for negative and positive 
values of v, respectively. Previous studies of low energy ion sputtering on Si surfaces 
suggest that the surface processes related to the second order derivative of the surface 
height are the curvature dependent erosion (CDE) rate |7J [T7J and ballistic smoothing 
[T8| [19j [20j [21]. The third term stands for a combination of relaxation processes while 
only ion-enhanced viscous flow (IVF) [6] are relevant for low energy ion sputtering on Si 
surfaces at room temperature [2TJ[22]. Equation [TJis linear and can be solved by Fourier 
Transformation, yielding for the (angular averaged) PSD function [201 [9] 

P(k,t) = P(M)exp(2r fc t) 
= — vk 2 — Dk 

where k denotes the spatial frequency, and the corresponding growth rate. In the 
case of v < 0, has a positive maximum at k c = ^—u/2D, leading to a correlated 
pattern with the characteristic frequency k c . 



2 ,,4 ( 2 ) 
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However, in contrast to the present results, the surface processes described above 
cannot directly be related to the presence of impurities at the surface. Shenoy et 
al [H] have shown that differences in sputtering yields and in surface diffusivities of 
different surface components can lead to modulations in topography and composition. 
In connection with the present experiments, a compositional inhomogeneity is indicated 
by the cross-sectional high-resolution transmission electron microscopy image of the 
dot patterned Si (100) surface shown in figure [2j The nanodots are crystalline bumps 
with an amorphous top layer of ~ 4 nm. The height and the periodicity of the dots 
are ~ 5 nm and ~ 50 nm, respectively. The intensity contrast is attributed to an 
accumulation of the metal impurities on top of the bumps. For a quantitative estimation, 
the darkened volume of ~2.5 nm height contains ~1.25xl0 16 Si atoms cm -2 . The areal 
density of metals concentrated on the crests can be estimated to be 2-3 times larger 
than the average total metal areal density of about 2xl0 15 cm~ 2 (see table 1). Thus, 
the composition is not inconsistent with that of a stoichiometric metal disilicide. This 
would confirm the conjecture by Ozaydin et al [H] that Mo silicides might be present 
during their ion sputtering experiments on Si surfaces. For the present system, a local 
decrease of the sputtering yield at the contaminated areas may be anticipated due 
to the larger surface binding energy of the silicide compared to pure silicon. Binary 
collision computer simulations using the TRIDYN program [23| 124"] confirm an decrease 
of the Si sputtering yield from FeSi2 vs. Si by ~ 10 %. Thus, once the non-uniform 
contamination is established, the differences of the sputter yield will indeed promote 
the pattern formation. There might be additional collisional effects which influence the 
pattern formation under metal contamination. The transition from holes to dots found 
by Sanchez-Garcia et al [13] was qualitatively attributed to a change in the shape of 
collision cascade due to the metal impurities in the near-surface region. 
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Figure 2. Cross-sectional HRTEM image of the correlated nanodots generated with 
metal incorporation. An amorphous top layer is seen which is darkened on top of the 
crests due to atomic number contrast. 

As an alternative or additional possible mechanism, local stress might influence 
the dot pattern formation during ion sputtering. Previous studies show that the 
stress generated in the near-surface region is large enough to influence the morphology 
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evolution [251 I2E1 EH EH] • Recently, Ozaydin et al [25] found a tensile stress which they 
attributed to Mo incorporation during ion sputtering. Such tensile stress could cause 
the growth of surface protrusions as a consequence of relaxation [29]. This might even be 
promoted by the non-uniform contamination described above. Nevertheless, it remains 
speculative if a stress mechanism could also contribute to the present findings. 

Thus, there is no conclusive information on the role of the metal incorporation 
neither from previous studies nor the present one. On the basis of the results presented 
above, we conclude that metal incorporation induces a new surface instability in addition 
to the curvature dependent erosion rate. Assuming that this new instability can be 
expressed as well by a second order derivative term as z/ m V 2 /i, we rewrite the growth 
coefficient r k = — (z/ CDE + v m + v^)k 2 — Dk 4 , where z/ CDE and 24, represent the curvature 
dependent erosion rate and ballistic smoothing, respectively. Correlated patterns are 
generated when isqde + v m + v^, < which means that the combination of the two 
instabilities dominate the fc 2 -dependent surface process and compete with fc 4 -dependent 
smoothing mechanisms. In the absence of metal incorporation no correlated pattern is 
formed on the surface, indicating that the instability induced by curvature dependent 
erosion rate is overwhelmed by the relaxation processes. It is worth to note that the 
PSD function of sample B shown in figure [1] exhibits a transition from a k~ 2 to a k~ A 
behavior at the spatial frequency of approximately 0.05 nm" 1 , as indicated by the two 
gray lines. This observation is in agreement with the asymptotic PSD function in the 
case of surface smoothing (r^ < for the entire spectrum) [201 E] 



where the PSD function is approximately proportional to A; 2 at small frequencies and 
k~ A at large frequencies. 

3.2. Smoothing of the pre-patterned surface 

For further investigation of the smoothing mechanisms, patterned samples were 
produced by ion sputtering with metal contamination as described above, and 
subsequently irradiated without metal contamination. 

Figure [3] show the AFM images from the original surface (a), which is similar to 
inset A in Fig. 1, and after post-sputtering with three different fluences of 1.0 x 10 17 
cm -2 , 2.0 x 10 17 cm" 2 and 3.0 x 10 17 cm" 2 ((b)-(d)), respectively. The evolution of the 
root mean square (rms) surface roughness obtained from lxl /im 2 AFM images and 
of the Fe surface content obtained from RBS are shown in figure H] (a). The morphology 
evolution can be divided into two regimes: At fluences smaller than 1.0 x 10 17 cm" 2 , 
the surface morphology does not visibly change, as shown by figure |3] (a) and (b). The 
roughness remains almost constant and even slightly increases. At fluences larger than 
1.0 x 10 17 cm" 2 , the nanodot pattern gradually decays with the inter-dot distance being 
unchanged until the surface becomes flat. An exponential decay of the roughness is 
observed. The smoothing effect in this regime is also confirmed by figure [3] (b)-(d). 




(3) 
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(a) (b) 




(c) (d) 




Figure 3. 1/imxl/im AFM images of Si(100) surfaces pre-patterned by 1200-eV Ar + 
ions at normal incidence with a stainless steel mask (a) and post sputtered with a Si 
covered mask at fluences of 1.0 x 10 17 cm -2 (b), 2.0 x 10 17 cm" 2 (c) and 3.0 x 10 17 
cm" 2 (d). 

The non-monotonic behavior in morphology evolution is attributed to the surface 
metal content. Based on the observation that the surface metal content decreases with 
ion ffuence (see Fig. 4 (a) for the Fe surface content), we assume that at the beginning of 
post-sputtering the metal atoms inherited from the pre-patterning process still remain 
on the sample surface although further metal incorporation is now prevented. Therefore, 
the surface instability is still large enough to compete with the smoothing mechanisms. 
Since the surface morphology of the pre-patterned samples has reached the saturation 
regime, the pattern does not change in the early stage of post-sputtering at fluences 
which are small compared to the pre-irradiation fluence. With post-sputtering being 
continued, the incorporated metal atoms are gradually removed from the surface, so 
that finally the smoothing mechanisms dominate the morphology evolution. In this 
sense, a threshold value of the metal content is expected for the roughening/smoothing 
transition. When the metal content is larger than this threshold value, the surface 
instability is large enough to trigger dot pattern formation. The present experiments 
indicate that this threshold value for a metal impurity of mass ~ 55 amu in Si is between 
0.6 x 10 15 cm -2 and 1 x 10 15 cm -2 , for Ar ion irradiation under the present conditions. 

Figure S] (b) shows the evolution of the PSD function in the decay regime. At 
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Figure 4. Post-sputtering of Si(100) surfaces pre-patterned by ion sputtering under 
metal incorporation, (a) Root mean square (rms) surface roughness W (squares) and 
surface content of Fe (circles) versus fluence $. The solid and dashed lines are plotted 
to guide the eyes, (b) PSD functions (calculated from 1 /imxl /im AFM measurement) 
of the samples post-sputtered with the indicated fluences. 



the first sight, we find smoothing over the entire spectrum with a rapid decrease of 
high frequencies and a successive decrease of lower frequencies. Finally, the morphology 
reaches a steady state. As expected, the steady-state PSD function shows a similar 
behavior as the one obtained after sputtering of a virgin sample without metal 
contamination (see Fig. 1, curve B). According to the discussion in section 13.11 the 
transient from a k~ 2 to a k~ 4 behavior in the PSD function corresponds to the asymptotic 
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solution (equation ([3])) in the case of surface smoothing. This is a strong indication that 
the relaxation processes include both k 2 - and /c 4 -dependent smoothing mechanisms. 

Surface smoothing during ion sputtering has been investigated in several studies. 
Vauth and Mayr [21] suggest that ion-enhanced viscous flow or ballistic smoothing 
dominates for high or low spatial frequencies, respectively, during keV ion smoothing 
of amorphous surfaces. Although Frost et al [20J observed only a /c 2 -dependent surface 
smoothing during ion sputtering on Si, they also found first evidence for a transition 
from a ballistic smoothing (oc k 2 ) to viscous flow smoothing (oc k A ) in the PSD function 
composed from two AFM images with different scanning size. IVF has been suggested 
in several studies as the dominant fc 4 -dependent smoothing mechanism during low- 
energy ion sputtering on Si surfaces [211 1221 |30] • The concept of ballistic smoothing was 
firstly put forward by Carter and Vishnyakov [18] as a result of the net displacement 
of the forward recoils moving preferentially along the ion beam direction. The ion- 
impact-induced downhill current mechanism proposed by Moseler et al [19] on the base 
of molecular dynamics computer simulation is considered phenomenologically virtually 
indistinguishable from Carter and Vishnyakov's proposal [31]. Zhou et al [32] observed 
fc 2 -dependent smoothing on AI2O3 surfaces during Ar + ion erosion at normal incidence. 
They attribute this effect to the ballistic atomic downhill current [19]. Accordingly and 
based on the discussion above, we suggest that ballistic smoothing and IVF are the 
two dominant relaxation processes during normal incidence Ar + ion bombardment of Si 
surfaces. 

4. Conclusions 

In conclusion, with 1200-eV Ar + ions bombarding Si (100) surfaces at normal incidence 
and below 100 °C, we have observed dot pattern formation and surface smoothing 
in the presence and absence of metal incorporation, respectively. Metal incorporation 
corresponding to less than one monolayer critically influences the roughening/smoothing 
behavior. In addition to the well-known curvature-dependent erosion rate, an additional 
surface instability is generated, which is consistent with laterally non-uniform sputtering 
due to preferential accumulation of metal atoms at the nanodots. In the case of 
smoothing, a transition from a k~ 2 to a k~ 4 behavior is observed in the asymptotic 
PSD function, which is ascribed to ballistic smoothing and ion-enhanced viscous flow 
as the dominant surface relaxation mechanisms. 
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